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We suggest a theory for a deformable and sliding charge density wave (CDW) in the Hall bar 
geometry for the quantum limit when the carriers in remnant small pockets are concentrated at 
lowest Landau levels (LL) forming a fractionally (v < 1) filled quantum Hall state. The gigantic 
polarizability of the CDW allows for a strong redistribution of electronic densities up to a complete 
charge segregation when all carriers occupy, with the maximum filling, a fraction v of the chain length 
- thus forming the integer quantum Hall state, while leaving the fraction (1 — z/) of the chain length 
unoccupied. The electric field in charged regions easily exceeds the pinning threshold of the CDW, 
then the depinning propagates into the nominally pinned central region via sharp domain walls. 
Resulting picture is that of compensated collective and normal pulsing counter-currents driven by 
the Hall voltage. This scenario is illustrated by numerical modeling for nonstationary distributions of 
the current and the electric field. This picture can interpret experiments in mesa-junctions showing 
depinning by the Hall voltage and the generation of voltage-controlled high frequency oscillations 
(Yu.I. Latyshev, P. Monceau, A.A. Sinchenko, et al, presented at ECRYS-2011, unpublished). 
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I. INTRODUCTION. 

Commonly, the sliding current of an incommensurate charge density wave (CDW) is driven by the electric field 
in the direction x of chains, which is originated by the applied voltage V and mediated by the normal current density 
jxn- New experiments [1] on mesa junctions in a high magnetic field (HMF) [U [3] show, for the first time, that the 
CDW can be driven also by the Hall voltage Vh originated by the normal current passing in a transverse direction, 
FigO The observations are quite counterintuitive with respect to the common notion. In a narrow channel cut 
transversely to the chains, the conductivity drops below the threshold rather than rising, Shapiro steps are quantized 
in voltage rather than in current, etc. 

Facing these challenges, we shall discuss here the Hall effect in the deformable CDW at the quantum limit of normal 
carriers. Our other studies show that the regime of classically strong magnetic fields (before installing the quantum 
Hall regime at lowest LLs) cannot reproduce the complex of experimental dependencies on B and sizes and the scale 
of the universal Josephson ratio for Shapiro steps. As a minimal model of the basic experimental geometry shown in 
Figj^ (left panel), we shall consider a Hall bar with the current and the magnetic field applied orthogonally in the 
interchain directions y,z. 

Small pockets of electrons can be present in the otherwise gapped CDW state. In a layered system at HMF the 
carriers are concentrated at one or a few lowest LLs forming fractionally (z/ < 1) filled QH state. As always, the 
Lorentz force pushes electrons to one edge x = 0 of the bar forming the slightly over/under charged boundary layers 
which give rise to the Hall voltage which strength Eh = VhI^x compensates for the Lorentz force. But now, the 
gigantic polarizability of the CDW reduces the electric field allowing for so strong redistributions of electronic densities 
which are unthinkable in usual circumstances. E.g. at low temperature TjlK, the experimental Vh ~ ImeV forces 
all carriers to completely fill a fraction u of the chain length, thus forming there the integer QH state, while leaving 
the fraction {1 — u) of the chain length unoccupied. 


II. OVERVIEW AND RESULTS. 

In a (quasi)-2D system at a HMF the electrons form essentially quantum states. But with transitions among the 
LLs being frozen, the kinetic energy leaves the game, then the density is determined classically by the electric potential 
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^(f). 

At edges of usual QH setups, the density distribution follows almost slavishly the strong Coulomb potential which 
has been build up already before the HMF was applied; it comes from pronounced finite size effects for the 2D 
electron gas embedded to the 3D geometry of the device, see im. But for electrons in pockets upon the CDW, 
deformations of the CDW phase are able to neutralize any charge so that electrons are allowed to form themselves 
highly inhomogeneous distributions. Here they are originated intrinsically, to provide the electric field opposing the 
Lorentz force. At low T, that extents to separation of the mean carriers’ density fii among a fully depleted - n = 0 
region where ^ and an overfilled one - n = n/v where ^ i.e. between the two, both nonconducting, 

regions where the same LL is either full or empty. 



FIG. 1: Left panel: The photo [T] of the working simplest mesa structure with comments. Dashed lines show the chains, 
arrows show the current which is forced to enter the long channel crossing the chains direction. Right panel: The scheme of the 
charge density segregation at T=0 into a region 0 < x < / of the maximal density in the IQH state and a fully depleted region 
I < X < Lx. Phase slips should take place at the boundary in the sliding regime. Jy is the bias transverse normal current. Jx 
and Jc are the compensating normal and collective currents in the chains direction. 

We consider a stack of plain rectangles of the size L^x Ly^ each composed of Ny chains with a spacing dy = Ny/Ly. 
Each chain segment bears a few A^i = fiiLx electrons. In the HMF they are supposed to occupy the lowest LL which 
can acquire as much as Nfl electrons. The mean LL filling is z/ = Ni/Nll = 27rA^ni/dy where A = ^J{hc/eB^ is the 
magnetic length. 

The profile for the density of carriers ni{x) is rii = /($ — ^F)ni/iy where f{S) is the Fermi distribution function at 
a given T and we have assumed that the LL energy hujd^^ is included to the definition of the Fermi potential energy 
The level is determined at the fixed number of particles by the condition 

J /(i>(x) - <!>F)dx = l = vF 

The in-plane current density in the direction y is given by properties of the Landau eigenfunctions in crossed electric 
and magnetic fields: 

3y{^) = “ ^f) 

Integration over x gives the overall y current as a function of only values of the boundary potentials: 

j exp[($F - ^(0))/r] + 1 

^ h exp[(i>^. - ^{L))/T] +1 ^ ^ 

In the limit T ^ 0 and if all states are filled (z^ = 1 which requires that both then we recover 

the standard result of the IQHE Jy = (^(0) — ^{L))e/h. In the limit T 0 for the partial filling u < 1, the states 
are all filled over a segment (0,/ = and all are depleted beyond, i.e. ^(0) < ^(0 = ^5 then 

Jy = (^(0) — ^i7)e/h depends explicitly only on the potential at one boundary. But this is only a part of the total 
voltage: the electric potential keeps growing also in the depleted region I < x < L^. The solution of the elementary 
electrostatic problem yields the value of ^{L) at the depleted boundary. Einally the mean current density reproduces 
the standard result for the classical Hall conductance r\j IIB: 

Jy = cr^cyi^iO) - ^{L)) , CTxy = ec^ 
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We have got that, in spite of the IQH filling of the populated area, because the division line between the areas is 
displaced with changing B. 

By now, there was no voltage drop in y direction which can be only the dissipative Vyy = pyyJy with a negligible 


Pyy — ^xx! {p XX'-’yy'\ ^ xy ) 

to the transverse one a^y^ In conventional QH devices it conies from localization of carriers (as the consequence of 
absence of their kinetic energy) at irregular potential of dopants. For electrons in NbSe3, with its \D 2D 3D 


jcjl-y. For that, the longitudinal conductivities Gxx , (^yy should be small with respect 


^ xx'-'yy 


< cr; 


xy 


for 


hierarchy, at least the interchain conductivity dyy is very small which ensures the condition a, 

domination of the transverse {-LB) currents j and of the Lorentz force. 

Allowing the CDW to start sliding brings about the normal counter-current = —Jc- The Lorentz force component 

appears in y direction giving rise to the additional voltage Vy. Its magnitude is given by the universal law similar to 

(1) with interchanged indices x and y (and the opposite sign). For the current per chain we get (the second line in 

(2) is written for T = 0): 


3 x 1 = dy^{x, y)f{£ {x, y))dy 

^ ($(x, 0) - L,)) = (2) 


Here AVy{x) is the mean interchain voltage drop; literarily, it is concentrated only over the filled stripe 0 < x < I 
and it is zero elsewhere while actually the source and the drain edges are expected to be equipotential. Following the 
experience [7] from the traditional IQHE, that can be corrected by a consistent 2D solution for the distribution of 
the carries density and the electric field. By now, we have imposed the perturbation from the x currents upon the 
already solved case of the bias y current. 

The last step is to remind that the CDW countercurrent, per chain, is jc = —edtpl^K and it requires for phase slips 
to annihilate with the normal current jxi at x = 0,/. Let each phase slip absorbs/releases M electrons per chain, 
then the repetition frequency is 


•’ 2wM 2Me ^ 2Mh ’ f e 

We arrive at the Josephson relation with an additional factor M; according to experiments, with a good accuracy 
M « 10 which provocatively coincides with the number of electrons per chain. It can happen that each phase slip 
swallows all electrons at the enriched end to deliver them at the depleted end. Such a strong event cannot be localized 
only near the ends and the whole dynamical problem needs to be modeled. 


III. NBSE3 IN HMF. 

Recall some relevant properties of the CDW in NbSe3. The most important feature is the gigantic dielectric 
constant e which source is the CDW polarization. There are several sources in the literature: e = 5 x 10^ at T=50K 
[8]; e = 2 X 10^ at T=42K, and e = 1 x 10^ at T=20K [9]; e = 5 x 10^ at T=22K [10]. We can extrapolate that to lower 
T ^ IK with a common wisdom (while not verified at low T) that e ^ 1/L^t which itself is known to be ^ exp(T/To) 
with To ~15.2K [TT1[T2|. Taking this law, or just points from plots in na, we expect a 20 times reduction from A2K 
to IK and 3 times from 22K] then we get the low T saturation at e ^ 10^“^. Notice that e ^ 10^ corresponds to the 
pinning length Ly « Ipk, then there will be no bulk pinning at all for our length scales. 

Now we summarize the facts related to normal carriers in NbSe3 at HMF. Hereafter all numerical estimations will 
be done for B = 20T. 

1. Zeeman splitting is Az ~ 2.32meV « 23K > Vh^ but still small to fully polarize the pocket of electrons which 
partial Fermi energy is 6meH, then 2/3 of electrons are available to form spinless pairs for an exchange with the CDW 
condensate. 

2. Magnetic length A ~ 5.64nm is very short in comparison with the chain length = 1pm. 

3. There are two types a of carriers: extrinsic a = e which are electrons and intrinsic a = i which are holes, with 
effective masses in the x^y plane, ujca = ^/(^a^^) are the cyclotron frequencies where huJce ~ ll.SmeH ~ 12K 
for electrons and an expected bigger value for lighter holes which are for sure in the truly quantum regime, since B^IT 
m- We shall suppose that bandwidths in the heaviest 2 : direction are smaller than huOc which statement is marginal 
for electrons and may be ensured for holes. Then electrons or at least holes form a 2D state within every x, y plane. 

4. The filling factor u of one LL is u = 7.8. Then only one-two LLs are necessary to accommodate the present ~ 10 
electrons per chain segment and only a small fraction of one LL is sufficient for expected ^ 1 holes. Precisely for 
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these numbers: for electrons, the n = 0 LL will be full and the n = 1 LL will be filled with a fraction ^ 0.3; for 
holes, the n = 0 LL will be filled with the fraction ~ 1/8. 


IV. THE MODELING. 

A. Conjections. 

We shall proceed with a minimalistic model which concentrates upon only one type of carriers at a zero LL; the spin 
splitting is not taken into account. We shall ignore the special low T states: the fractional QHE at special rational 
values of and the Wigner crystallization (its source, the long range Coulomb repulsion, is strongly suppressed in 
presence of the CDW). Electronic interaction will be taken into account only via the self-consistent slowly varying 
Coulomb potential 4>. This Hartry approximation was used effectively nn in theory of inhomogeneous IQHE near 
gates boundaries in semiconductors. 

We shall use the parameters as: chain length = 1/i/c, pinning length Ip = 0.1 fik (corresponding to e ^ 10^“^), 
temperature T = 0,0.1, 0.4 meV, magnetic field B = 20T, number of particles Ni = 1 per Lx (like for holes), filling 
factor of the LL as u = 0.1 and u = 0.3. We tested Vh = O.bmeV and Vh = 3meV corresponding to the threshold 
voltage and to the maximal one in experiments. 


B. Static regime of pinned CDW. 

Eor static distributions, the analytic treatment is possible at T=0 which may be applicable at the lower experimental 
limit T = IK. But for higher T > Vh^ we should use numerical solutions. The T = 0 solution confirms the full charge 
segregation which sets in when 

VH>Vg'= 2n{l-u)^l^N, 

edz dy 

The critical value ~ ImeV at e = 10^ is the energy of Coulomb interaction between the over- and under charged 
areas. 

Plots in Eigsj^show the numerically calculated distributions. 



FIG. 2: Static profiles for occupation numbers (cyan), currents (blue), electric potentials (black) and the CDW stress fic- The 
parameters are Vh = 3, Lp = 0.1; left panel: T = 0.1, u = 0.3; right panel: T = 0.4, u = 0.1 

• Cyan lines with fillings give distributions of occupation numbers - the Eermi function /((4> — 4>i7)/T). At Vh = 3mV 
the full depletion appears at the high-voltage end 0.8 < x < 1 when T = 1.1 A is low (cf. at T = 0 the depletion 
H = 0 would spread over 0.3 < x < 1); the IQH limit / = 1 is reached only near x = 0. At the elevated T = A.AK 
and smaller v = 0.1 the population is spread over the whole length, while it is strongly suppressed from / = 0.4 at 
X = 0 to f = 0.02 ai X = 1. 

• Blue lines give distributions of jy. They follow the rise of both and / towards x = 0. But the current can rise 
also towards depleted edge x = 1 where the high field can overcome the drop of the density. 

• Black lines give the monotonic potential 4>, with eE = — 4>' rising near both ends and suppressed, together with the 
current, in the middle. 










5 


• The brown line gives the CDW stress which is also the electrochemical potential jic of the condensate. The shape of 
jic is similar to that of $ but with a discrepancy A/i = ^ — which changes sign near the middle. That means the 
two reservoirs of electrons are not stable with respect to exchange of particles among them, provided the phase slips 
are allowed [Hill]. Then the CDW will give up its periods near the depletion region where A/i < 0 and will absorb 
pairs near the overpopulated end where A/i > 0. The discrepancy A/i is largest at the left end as we see it clearly at 
the plot for zi = 0.1, Vf/ = 3 where A/i reaches ^.2bmeV compatible with expectations for the phase slip voltage. 

Numerical results were obtained from the following equations. The one for the phase is the balance of the stress 
gradient and the pinning force Fp\ 




Me = ^ 


ttNf 


? — ^ 

^ ttNfLI 


(3) 


The potential satisfies the Poisson eq. with the total charge density per chain rifot = {q/tt n — n) 


+ m [f{<^>/T)/p -l]+q/n = 0, 

where Np = 2/iphvF), = 8e‘^/{hvFdydz), vq < Inm. Here q = ip' is the CDW wave number shift, —n is the 
background compensating charge density. 

With so small ro, eq. 0 is reduced to the constraint of the local electroneutrality 

Fm[f(^/T)/iy-l]Fq = 0 (4) 


Eqs. (^) were solved numerically with boundary conditions Here the first condition ensures the total electroneutrality 
while the second one fixes the Hall voltage. 


C. Time-dependent regime with depinning. 

We extend the above modeling in two respects. 

1. We shall model the dependence Fp{ip) in such a way that the linear law ^ sharply drops to zero dit cp > epp - the 
threshold phase displacement epp of depinning: 


Fp{(p) 


ip/jTTNFLl) 
lFex.p{k{(p - (pp)) 


( 5 ) 


This picture corresponds to a conventional first threshold field which looks to be always valid for NbSe3, while a more 
complicated two-step depinning is observed at low T in gapful CDWs, see [16]. 

2. We shall look for the time dependent evolution starting from the initial moment when the current Jy is turned 
on. We shall not limit the CDW current, given by the phase velocity u = dtp^ to vanish at the boundary. We shall 
suppose that the phase slips concentrated in narrow boundary layers will do this job. In the future, it should be 
possible to incorporate that to the model. 

Now, the eq. 0 is generalized as 


dxl^c = = ^ + dxip/{nNF) 


(6) 


Eqs. 0 , together with Q where q = dxP^ were solved numerically with the paired boundary condition for the electric 
field and separate conditions for potentials. 9a,4>(0,t) = dx^{Lx^t), ^(0,^) ~ = Vbtanh(t/r) < 0, ^{Lx) — ^f = 

Vi tanh(t/r) > 0 where the first condition ensures the total electroneutrality while others turn on the boundary voltage 
values gradually, over a short time r. Actually Vq and Vi are not fixed separately; rather the given current determines 
the relation Q among them. Notice that for low T when both | Vbl, ki ^ T it fixes only 4>(0, t) — ^f = Jyh/e leaving 
4>(I/a,,t) to be evaluated. 

Einally, we exclude 4> from eqs. to arrive at the eq. for the phase alone: 




tliNfT 


(ni - dxipli^) {dx<plT^ + ni(l - v)) 

^ILl 


+ 1 


(exp(fc((^ - ipp)) + 1) 


= idtp 


with Vi — Vq = Vh‘ Remarkably, such a nonlinear nonlocal partial-differential equation yields reliably a complicated 
spacio-temporal behavior described below. 
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FIG. 3: Space-time distributions of the electric field (left panel) and the phase velocity (right panel). 


The results presented in FigJ^were obtained for Lp = 0.1 , (pp = 0.1, k = 1000, = “I? = 2. Also we put 

7 = 1 which fixes the time scale. 

The left panel of FigJ^shows a plot of the electric field E{x^ t). It demonstrates two domain walls of the propagating 
depinning which collide and annihilate giving rise to a stationary sliding depinned state with essential but smooth 
variations over the length. The plot at the right panel shows that the phase velocity is not as sensitive to the passing 
walls but gives a strong, very narrow in time peak at their collision. 

The interpretation is that electric fields (created by electrons’ redistributions near the boundaries) exceed the 
threshold field while the bulk is still below that. The depinning penetrates into the pinned region by propagating 
sharp fronts of the electric field. An overshoot will be noticed if we plot also the distribution n(x, t)\ behind the fronts 
the occupation numbers are flat at 1 (the transient IQHE state) or zero (the full depletion). But later on the smooth 
profiles are established. 

By now, we studied only the single run (after switching on of the boundary potentials) assuming a source and a drain 
for the CDW current by phase slips. The implementation to the Hall bar and particularly to observed oscillations is 
anticipated to be the following, recall Sec.2. The CDW sliding is arrested at the boundaries by conversion of particles 
between the condensate and the normal states which job is made by phase slips. The rate is governed by difference of 
chemical potentials of the two reservoirs (notice the discrepancy between black and brown lines in plots of the Figj^. 
We expect that each impulse of the phase slip will reduce the conditions to the bare ones and the development similar 
to the above modeling will repeat itself sequentially. 

It seems that emergence of depinning waves in our modeling may be of an significant insight to other, more 
conventional experimental situations on sliding CDWs. 


V. OUTLOOK AND CONCLUSIONS 

In the bulk of sliding CDWs, the normal Jn and the collective Jc currents flow in the same direction while regulating 
their participation ratio depending on the geometry and the background m- A curious question may be posed: if 
the partial currents can flow to opposite directions? Even more, if they can be compensated to a limit when the 
total current is zero? It happens to be possible indeed in new circumstances when the on-chain currents are driven 
by the Hall voltage in a geometry with the closed circuit in a narrow x direction of the Hall bar, as described in this 
article and most probably already seen in the experiment [T] . Actually, the old experiments [18] on sliding under the 
temperature gradient alone should also be understood as the case of compensating currents. 

The picture of counter-currents (the collective and the normal ones) closing the loop by annihilating at the Hall 
bar boundaries looks rather fantastic and even suspicious as a perpetuum mobile. Still it is feasible, as well forgotten 
fundamentals of the conventional Hall effect in semiconductors m can tell us. In this old science there is a silently 
supposed and well forgotten feature of processes with more than one kind of carriers: presence of partial Hall currents 
for electrons and holes; the currents are compensated so that only the total one is zero. There, the currents of electrons 
and holes should be compared with the normal and the collective currents of the CDW. In both cases the equilibration 
among the carriers is hindered: by necessity of the energy relaxation across the gap in semiconductors or of making 
phase slips in the CDW. 

The CDW has a drastic influence already in the stationary regime, before the onset of sliding. Recall the theory 01^ 
[6] of inhomogeneous IQHE near gates boundaries in semiconductors. In those cases the stripe structure of electronic 
density appeared in response to the built-in electrostatic field spread from the gate boundary which is present already 









7 


at 5 = 0. In our case, the inhomogeneous distribution of concentration of electrons is caused internally by the Hall 
voltage. Any spread density profile deviating from the background charge of donors and the gate would be prohibited 
in a 2D electron gas in semiconductors because of the very high cost of Coulomb energy. In our case, the deformations 
of the CDW background allows to reduce the field by the gigantic dielectric constant e. The redistribution of charge 
caused by the Hall voltage alone, up to segregation into the overcharged fully filed IQH segment and the undercharged 
fully depleted one, is unthinkable in usual circumstance of the QHE. 
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